Phase diagrams of one-dimensional Bose-Fermi mixtures of ultra-cold atoms 
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We study the quantum phase diagrams of Bose-Fermi mixtures of ultracold atoms confined to 
one dimension in an optical lattice. For systems with incommensurate densities, various quantum 
phases, e.g. charge/spin density waves, pairing, phase separation, and the Wigner crystal, are found 
to be dominant in different parameter regimes within a bosonization approach. The structure of 
the phase diagram leads us to propose that the system is best understood as a Luttinger liquid of 
polarons (i.e. atoms of one species surrounded by screening clouds of the other species). Special 
fillings, half-filling for fermions and unit filling for bosons, and the resulting gapped phases are also 
discussed, as well as the properties of the polarons and the experimental realization of these phases. 
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I. INTRODUCTION 

In recent years the experimental control of trapped 
ultracold atoms has evolved to a level that allows to 
probe sophisticated many-body effects, in particular the 
regime of strong correlations. One of the most promi- 
nent achievements, the realization of the superfluid to 
Mott-insulator transition of bosonic atoms in an opti- 
cal lattice P, 0, triggered both experimental and the- 
oretical research into quantum phase transitions of ul- 
tracold atoms in optical lattices. The advances in trap- 
ping and cooling techniques, as well as the manipulation 
of atom-atom interactions by Feshbach resonances Q, 
created the possibility to study many-body systems in 
the quantum degenerate regime in a widely tunable en- 
vironment. Of particular interest are those systems that 
resemble solid-state systems, for example, mixtures of 
ultra-cold bosonic and fermionic atoms Q, which have 
recently become accessible through the development of 
sympathetic cooling |5]. On the theoretical side several 
phenomena have been proposed in the literature like pair- 
ing of fermions Q , formation of composite particles Q , 
spontaneous breaking of translational symmetry in opti- 
cal lattices |{|. However, the approach used in most of 
these studies, which relied on integrating out the bosonic 
degrees of freedom and then using a mean-field approach 
to investigate many-body states [6J, becomes unreliable 
in the regime of strong interactions, in particular, in low- 
dimensional systems. A non-perturbative and beyond- 
mean- field investigation is therefore of high interest. 

In onc-dimcnsional (ID) systems, on the other hand, 
many physical problems can be solved exactly by various 
mathematical methods because of the restricted phase 
space and higher order geometric symmetry. One of the 
most successful examples is the Luttinger liquid model 
for ID electron systems 11 CM . The pioneerin g w ork of 
Tomonaga [ll[ , Luttinger jl2j and Haldane [H, 14 1 along 
with many others has produced an essentially complete 
understanding of the low energy physics of these sys- 
tems (within the so called bosonization approach). Re- 
cently these theoretical methods were also applied to sys- 



tems of ultra-cold atoms [HI, [ljl [TtJ , where the bosonic 
or fermionic atoms are confined in a highly enlongated 
magnetic-optical trap potential and approach an ideal 
ID or quasi-lD system. Such ID or quasi-lD enlongated 
trap potentials can be easily prepared either in a tra- 
ditional magnetic-optical trap [lj| , magnetic waveguides 
on microchips [l!|, or in an anisotropic optical lattice 
[20| . However, since there is no true long-range order in 
a homogeneous ID system in the thermodynamic limit, 
the one-dimensional "quantum phases" we will refer to in 
this paper are actually understood in the sense of quasi- 
long-range order (QLRO), i.e. the correlation function of 
a given order parameter (0(x)) has an algebraic decay at 
zero temperature, (O(x)O(O)) ~ |x| _2+Q , where (■ • •) is 
the ground state average, x is the distance, and a is the 
scaling exponent associated with that order parameter. 
As a result, the dominant quantum phase is determined 
by the largest scaling exponent a > (i.e. the slowest 
algebraic decay of the correlation functions of a given or- 
der parameter). At the phase boundaries that appear in 
the phase diagrams in this paper, the scaling exponent 
of the corresponding type of order becomes positive, i.e. 
the corresponding susceptibility switches from finite to 
divergent. 

In this paper we extend and elaborate on our previous 
work [TU on one-dimensional (ID) Bose-Fermi mixtures. 
We give a detailed derivation of the effective low-energy 
Hamiltonian of a ID BFM within a bosonization ap- 
proach [3, EH . After diagonalizing this Hamiltonian we 
calculate the long-distance behavior of the single particle 
correlation functions. We introduce a fermionic polaron 
(/-P) operator, constructed out of a fermionic operator 
with a screening cloud of bosonic atoms, and determine 
the single-particle correlation function, as well as the cor- 
relation functions of various order parameters, which are 
used to construct the phase diagram. For a BFM of one 
species of fermions with one species of bosons, we find 
a charge density wave (CDW) phase that competes with 
the pairing phase of fermionic polarons (/-PP), giving a 
complete description of the system as a Luttinger liq- 
uid of polarons. When the boson-fermion interaction 
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strength is stronger than a critical value, we observe 
the Wentzel-Bardeen instability (WB) region, where the 
BFM undergoes phase separation (PS) for repulsive in- 
teraction or collapse (CL) for attractive interaction. 

We also study several special cases when the density 
of the fermions or the bosons is commensurate with the 
lattice period. The low energy effective Hamiltonian of 
these cases cannot be diagonalized but instead has to be 
studied by using a renormalization group (RG) method. 
We obtain a charge gapped phase for half-filling of the 
fermions and for unit filling of the bosons, and study how 
these phases are affected by the presence of the other 
species. 

We also apply the bosonization method in an anal- 
ogous approach to spinful fermions mixed with a sin- 
gle species of bosonic atoms. We find a rich quan- 
tum phase diagram, including spin/charge density waves 
(SDW/CDW), triplet/singlet fermionic polaron pairing, 
and a regime showing the Wentzel-Bardeen instability. 
The similarity with the known phase diagram of ID in- 
teracting electronic systems [l(|, again suggests that a 
f D BFM is best understood as a Luttinger liquid of po- 
larons. We also discuss the close relationship between 
the polarons constructed in the bosonization approach 
in this paper and the canonical polaron transformation 
usually applied in solid state physics. 

This paper is organized as follows: In Section [TTJ we 
present the microscopic description of a BFM and its 
bosonized representation. In Section Ifffi we calculate 
the scaling exponents of single particle correlation func- 
tions as well as the correlation functions of various or- 
der parameters. Quantum phase diagrams obtained by 
comparing these different scaling exponents arc presented 
in the anti-adiabatic regime, i.e. the phonon velocity is 
much larger than the Fermi velocity. In Section IIV1 we 
present the results of two commensurate filling regimes: 
half-filling of fermions, and unit filling of bosons. In Sec- 
tionfVl we study a BFM of spinful fermions with SU(2)- 
symmetry and determine the rich phase diagram of this 
system. In Section fVTl we discuss questions regarding the 
experimental realization of the quantum phase diagrams 
presented in this paper, and we summarize our results 
in Section IVIII In Appendix [5] we give a more technical 
derivation of the scaling exponents of various operators. 
In Appendix [Bl we discuss the notion of polarons used in 
this paper, and compare it to the those used in solid state 
systems. The effects of a local impurity in the system are 
discussed in Appendix [Q 



II. BOSONIZED HAMILTONIAN AND EXACT 
DIAGONALIZATION 

In this section we derive and discuss the effective 
low energy Hamiltonian of a ID BFM in the limit 
of large bosonic filling 3> i/f, where v b /f is the 
bosonic/fermionic filling fraction) and fast bosons (u& 3> 
Vf, where Vb is the bosonic phonon velocity, and Vf is 



the Fermi velocity). We also mention how this study 
can be extended to the regime of comparable velocities 
(v f ~ v b ). 

We consider a BFM in an anisotropic optical 
lattice, where the lattice strengths for bosons (b) 
and fermions (/) can be expressed as V&(/)(r) = 
J2 a =x.,y,z V b(f),aSm 2 (kir a ) with h = 2ir/\ being the 
laser wavevector and A being the wavevector of the laser 
fields. Throughout this paper we will use A/2 = 1 (which 
corresponds to the lattice constant of the optical lattice) 
as the natural length scale. In order to create an effec- 
tively ID lattice, we consider V&(/vii <C Vb(/vj_, where 
— Vb(f),x is the lattice strength along the longitu- 
dinal direction and V(,(/),X = = Vuf) >z is along the 
transverse direction, so that the single particle tunneling 
between each ID tube is strongly surpressed. We note 
that independent tuning of the optical lattice strength 
for these two species of atoms can be achieved even when 
only a single pair of lasers provides the standing beam in 
each direction. This can be done due to the fact that 
the strength of an optical lattice for a given atom is pro- 
portional to Q 2 R /(v — AE) [2l|, where flu is the Rabi 
frequency (which is proportional to the laser intensity), 
v is the laser frequency, and AE is the resonance en- 
ergy of that atom — therefore for two different species 
of atoms (different AE), their effective lattice strengths 
can be tuned independently over a wide range, by simul- 
taneously changing the laser intensities via £Ir and the 
laser frequency v. 



A. Hamiltonian 

For sufficiently strong optical potentials the micro- 
scopic Hamiltonian is given by a single band Hubbard 
model [l|: 

H = - (tbb\bj + tfflf^j - Y (fJ>f n f,i + M6«6,*) 

(ij) i 

+ ~Y ^2 n t>A n bS - 1) + U b f Y n b,i n f,ii (!) 
i i 

where nb,i = b\b i and n/,, = f\f i are the boson and 
fermion density operators, and /ib/f are their chemical 
potentials. The tunneling amplitudes and the par- 
ticle interactions Ub and Ubf can be calculated from the 
lattice strengths and the s-wave scattering lengths: 

t b (f) = (4/V^)n 3 ( /* || ^ 6(/) exp(-2^ i|| ) (2) 
U b = (8/n) 1 / 2 (2ira b b/\)Eb(Vb,\\V 2 1 _) 1 / 4 (3) 
U bf = (8M 1 / 2 (2na b f/X)Ebf(V ba V h l^ 4 (4) 

We defined V r 6(/),||/± = V bU) ^/j_/ E b{f} and = 
(4H,||(±)V/ ) ||(X))/((n;^ ) + V f f {±) ) 2 ). E bif) = 
(2ir/ X) 2 /2nib(j) are the recoil energies of the two 
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atomic species, and Ebf is given by £y = (mj + 
rrif )/Am b mf (2n / A) 2 . 

For the convenience of the subsequent discussion we 
also give the Hamiltonian in momentum space: 



H = H<®+Hf> + H hJ 



(5) 



Here, refers to the bosonic sector without the cou- 
pling to the fermions, and is given by: 

Hi° } = E( £ M ^ I^Ah + p{ kPb , k (6) 



L is the ID system length, eb.k is given by tb,k = 
-2tbCOsk, and pb,k is defined as pb,k — H p b P+k b p . The 
free fermionic sector - without the coupling to the bosons 
- is given by: 



k 



(7) 



where ej^ k is given by e/^ = ~2tfCOsk. The interaction 
between these two sectors can be written as: 



k 



H bf - —f- 2^P'f,kPb,k 



(8) 



with p fik = Y.pft+kfp- 



B. Low energy effective Hamiltonian 

To determine the effective low energy Hamiltonian, we 
first consider the effective Hamiltonian of weakly inter- 
acting bosons, , without coupling to the fermions. 
Following the standard Bogoliubov transformation [2l|, 
Hll, we obtain the following effective low-energy Hamil- 
tonian for the bosonic sector: 



(9) 



k^O 



where the (3 k arc the Bogoliubov phonon operators and 
Lo b ,k = y/Jebjs - £6,o)(eb,fe - e&,o + ^U b Vb) is the phonon 
dispersion. The phonon operators 0k are related to 
the original boson operators by (3 k — u k b k — v k b_ k , 



where 



= 1 



(efc,fc — £fc,o + U b Vb + Wb,k)/2Ub,k 



and UkVk = —UbVb/^b,k are the coefficients of the Bo- 
goliubov transformation [2l[. We note that, although 
this tranformation is considered as a treatment of the 
quadratic fluctuations around a meanfield approxima- 
tion, it is still applicable to ID systems where no true 
condensate exists even at zero temperature (see [3^|V In 
the long wavelength limit, we have u>b t k — Vb\k\, where 
Vb = \J2tbUbVb is the phonon velocity. 

For the fermionic sector we proceed along the lines of 
the Luttingcr liquid formalism, established in solid state 

physics [m LO, El, - We linearize the noninteracting 



fermion band energy around the two Fermi points (at 
ikf, with kf being the Fermi wavevector, kf = TTVf), 
and split the fermion operator into a right (R) and a left 
(L) moving channel. One can show that such a linearized 
band Luttinger model [11] is effectively equivalent to a 
bilinear bosonic Hamiltonian in the low energy limit: 



H 



(0) 

/ 



(10) 



where Bk is a bosonic (density) operator defined as 
Bk = tVSEp/L+p/fl,, fOT k > and B k = 



*y \k\L fh,k+pi L,p 

right/left movers and 



for k < 0. fn/L,k are the 



2t/sm(fe/) 



(11) 



is the Fermi velocity. 

Finally we will address the interaction between the 
bosonic and the fermionic atoms. In the low-energy limit, 
there are two kinds of scattering to be considered. One 
is the scattering by exchanging small momentum in both 
the fermionic and the bosonic sector. The other one is 
the scattering by exchanging a momentum of 2k f to re- 
verse the direction of a fermion moving in ID. For the 
first term, only the long wavelength density fluctuations 
are important so the corresponding interaction can be 
expressed by the same transformation used above: 



H 



(i) 
bf 



Y,9\k\([3l+f3-k)(B k +Bl k ). 

k^O 



where g and Kb are given by 



K,, 



TX\j2t h VblVb 



(12) 



(13) 
(14) 



Kb is the Luttinger parameter for bosons as will be dis- 
cussed further in the next section. 

For the second case, low energy excitations occur in 
the backward scattering channel, and therefore cannot 
be included in Eq. (|12p . shown above. Such a backward 
scattering term, in general, leads to a non-linear and non- 
diagonalizable term in the Hamiltonian. We therefore 
consider a certain parameter regime in which the back- 
ward scattering term becomes tractable. In this and the 
next section . we will consider a BFM with large 
bosonic filling (i.e. i/j 3> ^/), and fast boson limit (i.e. 
Vb 3> Vf). Experimentally, this limit can be achieved by 
choosing tb S> t/, and Vb ~ 1 — 3, for typical atomic in- 
teractions. For simplicity, we will also assume that both 
the bosonic and the fermionic filling fraction are not com- 
mensurate to the lattice period or to each other. We will 
show that the Hamiltonian of such a system can be solved 
exactly even in the presence of backward scattering. 

In the limit of large bosonic filling and fast boson ve- 
locity, it is easy to see that the backward scattering of 
fermions is mainly provided by the forward scattering of 
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bosons (because 2kf = 2m/ f <C 2fcb = 27r^) , so that the 
effective interaction Hamiltonian in this channel can be 
written as 



H, 



(2) 
bf 



1 

7l 



52 fc/ (4+ME4 +fc /p( 15 ) 



\k\~2kf 



which is bilinear in two different bosonic operators and 
therefore can be diagonalized exactly. For the con- 
venience of later discussion and calculating the corre- 
lation functions, in the next section we will use Hal- 
dane's bosonization representation to express the effec- 
tive Hamiltonian shown above. 



where the coupling is given by g k — Ubf[vb(^b, k — 
£b ] o)/27ro;b ] / c ] 1 / 2 . The next step is to separate the 
non-interacting phonon field, Eq. ([9]), into a low en- 
ergy and high energy part, H = T,\k\<k f UkPkPk + 
Y2\k\>k ^kPt-Pk) and then integrate out the latter for 
|fc| ~ 2k f, together with the backward scattering term 
in Eq. (jT5]). Since we assume vj, 3> Vf, we can obtain 
an effective interaction between fermions within an insta- 
neous approximation. After some algebra, we can obtain 
such an effective fermion-fcrmion interaction to be: 



H 



(i) 



i \k\~2k f v q 



1 92k t 



L LU2k s 

1 9 2 2kf 



} [n R> k>n Ltk > + n L , k <n R , k '} (16) 



where n R/L>k = Y, p fn./L,p+ k fR/L, P is the density opera- 
tor for the right/left moving electrons. We therefore find 
that the original backward scattering obtained by inte- 
grating out the boson field now becomes a forward scat- 
tering term with repulsive interaction between the left 
and right movers. Therefore we can apply the bosoniza- 
tion approach for Luttinger liquids used in Eq. (fit))) and 
obtain 



H { p = Y,G\k\(B k B- k + BtBl k ) 



(17) 



fc~0 



C. Effective Hamiltonian in Haldane's bosonization 
representation 

We now bosonize both fermions and bosons on equal 
footing by introducing the phase and density fluctuation 
operators [14 , Il6j through the definitions: 



r -re"' (20) 



f{ x ) = [vf + Ylf] 112 « (2m,1)ie 

ra— — oo 

oo 

b(x) = K + n b ] 1/2 e 2rnl0b e 



(21) 



m— — oo 



where we use x to denote a continuous coordinate if there 
is no lattice potential in the longitudinal direction or a 
discrete coordinate in the presence of a lattice. H b /f(x) 
and & b /f(x) are respectively the density and phase fluc- 
tuations with Ob// = Tris b /fX+TT f dyHf,/f(y) accounting 
for the discreteness of atoms along the ID direction. The 
density and phase fluctuations satisfy the commutation 
relations: 



[lL b /f(x),$ b/ f(x')] =i8{x-x'). 



(22) 



where 



The advantage of using the density-phase fluctuation rep- 
resentation for boson-fermion mixture is that it treats the 
fermions and bosons in the same way. We do not have 
to use a Bogoliubov transformation for the bosonic field 
and a Luttinger liquid formalism for fermions separately 
as has been done above. As has been established in the 
literature (e.g. [H], Ell), the effective low energy Hamil- 
tonian of the bosonic and the fermionic sector can be 
written in terms of the fields ^l b /f and f as: 



G 



9zkf_ 



(18) 



Therefore the sum of the terms given in Eq. (jTTJJ) , 
([12"]) and (fT7|) is our final effective low energy Hamiltonian 
within the limit of large bosonic filling and fast boson 
velocity. The total Hamiltonian then becomes: 

H = J2^b, k (3t( 3 k + J2 v f\ k \ B l Bk 



k^O 



k^O 



Y,9\k\{Pl + P-k){B k + Bl k ) 



k^O 



■Y,G\k\(B k B_ k + BlBl h ), 



(19) 



(0) 



H, 





^ dx 


'Kb 

7T 




f dx 




2 j 


7T 



K, 



n 2 



K 



(23) 
(24) 



where we have Kt = 1 due to the absence of s-wave 
short-range interaction between fermionic atoms. Note 
that for intermediate values of the bosonic interaction 
Ub, the phonon velocity Vb and the Luttinger exponent 
Kb are renormalized by higher order terms. They can be 
obtained by solving the model exactly via Bethe-ansatz 
and are found to be very well approximated by 



Vb 



fe~0 



— V7 

ml v 



2tt 



1/2 



(25) 



5 



and 



K h 



V7 



1 



2tt 



-1/2 



(26) 



where is the effective boson mass in the presence of 
the lattice potential [1H and v s ~ ^ is the superfluid 
fraction. 7 = mlUb/v s is a dimensionless parameter, 
characterizing the interaction strength of the ID boson 
system. To leading order in the interaction (in the weakly 
interacting limit 7 — > 0), u& and Kb are the same as 
before. 

The density and phase fluctuations, defined in Eqs. 
([20]) and (f2Tj) , can be related to the bosonic operators f3 k 
and Bk by (see [HI ]) 



*i(x) 




n(k)e tkx ( 



B k + B 



-k 



(27) 
(28) 
(29) 



fc^O 



- S gn(k)e ikx (B fc -Bl fc )30) 



which can also be seen by comparing Eqs. (|23|) -([24 |) with 

Eqs. ©-cni) 

The long wavelength limit of the boson-fermion inter- 
action and the backward scattering Hamiltonian shown 
in Eqs. (fT2"|) and ([!"?]) can also be expressed in terms of 
the density and phase fluctuations operators as follows: 



H, 



H 



(i) 



= g / dx [, 2 U 2 



(31) 
(32) 



Therefore the total low-energy effective Hamiltonian is 
given by the sum of Eqs. ([23]), (ED), (31]), and (32]): 



(o) 



H, 



(o) 



(i) 
ft/ 



H 



(i) 



(33) 



which can be diagonalized easily as described below. 

Before solving the effective Hamiltonian, we note that 
similar problems have been investigated in the context 
of electron-phonon interactions in one-dimensional solid 
state physics in the literature, both with [25j, |27j or with- 
out [2y, [23I backward scattering of electrons. In the 
previous studies, the backward scattering of fermions 
was treated in bosonized form before integrating out the 
phonon field. The resulting effective interaction obtained 
after integrating out the high energy phonon field is of 
the form [25j] 



Hf ] cx J J d 2 rd 2 r'D {r - v 1 ) cos [26/(r)] cos [26/(r')] 



(34) 



where r = (x,w/r) is a space vector in 1+1 dimension 
and Do(r — r') is the phonon propagator. (Note that 
the symbol of density and phase fluctuations used in this 
paper and in Eq. ([M]) is different from Rcf. 25].) The 
instantaneous approximation of the phonon propagator 
is given by Dq(y — r') cx S(r— r'). However, Eq. (|3"4"|) does 
not seem to reduce to a bilinear term as shown Eq. (|32p . 
Such superficial disagreement can be resolved by taking 
the normal ordered form of cos (29) carefully, as treated 
by Sankar in Ref. [30]. In our treatment we keep the 
fermionic form after integrating out the phonon field and 
rearranging the order of fermions before using bosoniza- 
tion, which is a technically simpler and unambiguous pro- 
cedure. It shows that within the instantaneous limit the 
ID BFM system (and the related electron-phonon sys- 
tem) can be described by a bilinear bosonized Hamilto- 
nian, which then can be diagonalized exactly, even in the 
presence of backward scattering. 

We note that the generic form of the effective Hamil- 
tonian (|33[) extends beyond the parameter regime as- 
sumed in this section. If we drop the assumption of fast 
bosons (vb ^> Vf), and rather assume these velocities to 
be comparable, a renormalization argument of the type 
presented in [251 ] shows that the resulting Hamiltonian 
can also be written as (|3"3"|) . Then, however, the expres- 
sions for the effective parameters shown in this section 
would not hold anymore, and would need to be replaced 
by parameters resulting from the RG flow. 



D. Diagonalized Hamiltonian 

We now go on to diagonalize the full effective low en- 
ergy Hamiltonian (Eq. I|33p) by using the following linear 
transformation [261 ] : 

n fc = (5in A + <5 2 n a , $ 6 = e^ A + e 2 $ Q , 
Tl f = /3 1 n A + /? 2 n a , $/ =7i*A+72$a, (35) 

where Ha/ a and &A/a are the density and phase opera- 
tors of the two eigenmodes, and the coefficients are given 
by: 




(36) 



cos ip. 



Here we have e e = {{v f - 2G)/(v f + 2G)) 1 / 4 , = y/K~ b , 
and tan 2ip = Ag^bVf) 1 / 2 / (v 2 — v 2 ). The diagonalized 
Hamiltonian 

1 



j=a,A 



dx 



nIL j (x) 2 + -(d x $ j (x)y 



(37) 



has the eigenmode velocities, va and u a , given by 



V l/A 



1 



1 



6 



where Vf = (v 2 — 4G 2 ) 1 / 2 and g = ge e , e e is given by 
e° = ((vf-2G)/(vf + 2G))V\ 

From Eq. p8[) we note that when the fermion-phonon 
coupling g (proportional to Ubf) becomes sufficiently 
strong, the eigenmode velocity va becomes imaginary, 
indicating an instability of the system to phase separa- 
tion (PS) or collapse (CL), depending on the sign of Ubf 
[l7L [23l | . This is the so-called Wentzel-Bardeen instabil- 
ity, and has been studied before in ID electron-phonon 
systems [l7[ 



III. SCALING EXPONENTS AND QUANTUM 
PHASE DIAGRAMS 

In this section, we will calculate the scaling exponents 
of various quasi-long-range order parameters to construct 
the quantum phase diagram. The diagonalized form of 
the Hamiltonian, which we presented in Section [II Dl al- 
lows an exact calculation of all correlation functions (see 
Appendix |XJ). For the later discussion, we define the fol- 
lowing "scaling exponents" to mimic the ones in standard 
single component bosonic (l6j or fermionic systems [Io| : 

K fj = (i\ + pi K- 1 = 7l 2 + 72 2 (39) 
K S = Sj+Sl K~ 1 = 4 + 4 (40) 
Kps = P161+P2S2 K~l =7iei +7 2 e 2 . (41) 

Note that since v a > va according to Eq. (f3"5|) . we have 
Kjss < and K~} > 0. From Eq. (|55|l it is apparent that 
Kp and K~ x are the scaling exponent associated with the 
fcrmion density and phase fluctuations, respectively, and 
K$ and K~ x are associated with the boson density and 
phase fluctuations. Kps and are related to fluctu- 
ations of mixed components, and appear in the scaling 
exponents of products of fermionic and bosonic opera- 
tors. 



A. Single particle correlation functions and 
polaron operators 

In this section we consider the single-particle correla- 
tion functions of the system. For the bare bosonic and 
fermionic particles we find (b(x)b'(0)) ~ |a;| _ 2^e and 
(f(x)P(O)) ~ cos(fc / a;)|x| _ 3(^+^ 1 ), as discussed in 
Appendix [X] The scaling exponents that appear in these 
correlation functions have been renormalized by the cou- 
pling between the bosons and the fermions. 

Motivated by the polaronic effects in electron-phonon 
systems, we next consider a class of operators that seems 
to be most appropiate as the elementary operators of 
the system. In BFMs each atom will repel (attract) the 
atoms of the other species in its vicinity, due to their mu- 
tual interaction, resulting in a reduced (enhanced) den- 
sity around that atom. To describe particles dressed with 



screening clouds of the other species we introduce the fol- 
lowing composite operators (polarons): 

f x (x) = e- iX *»Wf(x), (42) 
b n (x) = e-">*f {x) b{x), (43) 

with A and r\ being real numbers, describing the size of 
the screening cloud. The correlation functions of these 
new fermionic and bosonic operators can be calculated 
to be (cp. Appendix lA|l 

<A(z)/](0)> ~ cos^M-^^ 1 - 2 ^^ 2 ^ 1 ) 

(44) 

(b v (x)ll(o)) ~ i^-i^W^- 1 -^- 1 ). (45) 

Treating A and rj as variational parameters, we observe 
that the exponents of the correlation functions are min- 
imized (i.e. the correlation functions have the slowest 
decay at long distances) by taking A c = K e /K ie and 
T] c = K 1 fK lli . To gain some insight into the quantity A c 
we consider the limit of weak boson-fermion interaction 
(Ubf — > 0) and obtain approximately A c ~ Ubf jUb (for 
further discussion, see Appendix |B|) . This result can be 
interpreted by a simple density counting argument as fol- 
lows: we imagine a single fermionic atom interacting with 
—A bosons in its vicinity. The potential energy of such a 
configuration can be estimated as —Ubf\+Ub\ 2 /2, which 
is minimized by taking A = Ubf/Ub- Since the polaronic 
cloud size is not a well-defined quantum number, A c and 
r\ c need not be integers and depend on the boson-fermion 
interaction strength continuously. 

Taking the optimal values of A c and r\ c in Eqs. ([42]) - 
(|43l and using the following two identities: K^ 1 = K~ x — 

K e K~ 2 and KJ 1 = K~ x — K 7 K~ 2 , we can show that 

(fxMflm - cos(k f x)\x\~^ +K ^ (46) 

(M^K Q )> ~ (47) 

In other words, the long-distance (low energy) behavior 
of both the fermionic-polaron (/-P, defined in Eq. ([42]) 
with A = A c ) and the bosonic-polaron (6-P, defined in Eq. 
(|43|) with 77 = rye) are characterized by a single scaling 
parameter, Kp and Kg, respectively, in contrast to the 
scaling behavior of the bare atoms. As we will see in 
the next section, the many body correlation functions 
of operators composed of such fermionic polarons also 
scale with the same single parameter, Kp, leading to the 
conclusion that the ID BFM system is best understood as 
a Luttinger liquid of polarons, instead of bare fermions. 

B. Correlation functions of order parameters 

It is well established that there is no true off-diagonal 
long range order in ID systems, due to strong fluctua- 
tions. Therefore, in this paper, the ground state is char- 
acterized by the order parameter that has the slowest 
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FIG. 1: (a) Scaling exponents of various order parameters 
as a function of longitudinal FP coupling g for a BFM with 
spinless fermions. Parameters are chosen to be: vt/vf — 3, 
Kb — 5 and G/vj = 0.1. Different curves correspond to 
the 2fc/-CDW, /-polaron pairing field (/-PP), bare fermion 
pairing field (BFP), 6-polaron operator (6-P), and bare boson 
operator (BB). Note that operators constructed with polarons 
(i.e. &-P and /-PP) always have larger exponents than their 
counterparts constructed with bare atoms, (b) Global phase 
diagram of the same system in terms of effective forward scat- 
tering between fermions, G and the fermion- phonon coupling 
g for Vb/vf =5 and Kb = 10. Three different phases, CDW, 
/-PP, and phase separation (PS) are dominant successively 
from weak to strong fermion-phonon coupling strength. 

long distance decay of the correlation function [To| , i.e. in 
the sense of QLRO. Transforming to momentum-energy 
space, this is equivalent to finding the most divergent sus- 
ceptibility in the low temperature limit: if, at T = 0, the 
correlation function of a certain order parameter 0(x) 
decays as (O(x)O^(O)) ~ l/\x\ 2 ~ a for large x, the fi- 
nite T susceptibility will diverge as x(T) ~ 1/T a . Here 
we defined the scaling exponent a, which is positive for 
a quasi long range ordered state (i.e. x(T) diverges for 
T — ► 0), and negative if no such order is present. Across a 
phase boundary the susceptibility will switch from finite 
to divergent. Since we are interested here in the limit 
of high phonon velocity (vb S> Vf), only the fermionic 
degrees of freedom need to be considered for the low- 
energy (long-distance) behavior of the correlation func- 
tions. The bosonic sector is always assumed to be in a 
quasi-condensate state. (We will relax these assumptions 
below.) 

Similar to a LL of spinless fermions, we first con- 
sider the order parameters of the CDW phase, defined as 
Ocdw — fi,( x )fR( x )t an d °f the p-wave pairing phase, 
given by Obfp = fh{x)fR{x), where the abbreviation 
BFP stands for bare fermion pairing. The scaling expo- 
nents for the CDW and the BFP phase can be calculated 
(see Appendix [A"f to be: 

a C DW = 2-2K p (48) 
a BFP = 2-2K- 1 (49) 

In Fig. [1] (a), we show the calculated scaling exponents 
for CDW and BFP phases as a function of the long- 
wavelength fermion-phonon coupling strength g, which 
is proportional to Ubf- The other parameters are given 
by Kf, = 5, Vb/vf = 3, and G/vf — 0.1. One can see 



that CDW ordering becomes dominant (chcdw > 0) for 
g/vf < 0.55, while bare fermion pairing is dominant 
{oibfp > 0) when g/vf > 0.65. In the intermediate 
regime (0.55 < g/vf < 0.65) none of these two phases 
have quasi long range order, and therefore one might con- 
clude that in this regime the system displays a metallic 
phase. A similar result has been discussed in the context 
of ID electron-phonon systems, [25l. 1271. l28j. 

However, as mentioned in our earlier publication [l5l |. 
the above analysis is incomplete. An indication for this 
is given by how a single impurity potential affects the 
transport of the ID BFM (see Appendix [Cj . Following 
the renormalization group analysis by Kane and Fisher 
[37l | . we show that a single weak impurity potential is 
relevant when olcdw = 2 — 2Kp > and becomes irrele- 
vant whenever the system is outside the CDW regime, i.e. 
when acDW = 2 — 2Kp < 0. This result is inconsistent 
with the previous CDW-metal-BFP scenario, because a 
metallic state would not be insensitive to an impurity 
potential in ID. More precisely, the fact that the impu- 
rity potential becomes irrelevant even outside the BFP 
phase indicates that there is another superfluid-like or- 
dering, which cannot be described by the pairing of bare 
fermions only. 

Above observation motivates a further search for the 
appropiate order parameters. Among the many types 
of operators that can be considered, we find that the 
operator 0/-pp = fL,\ c (x)fR,\ c (x), which describes a 
p-wave fermionic polaron pairing phase (/-PP), shows 
dominant scaling outside of the CDW regime. In Fig. [5] 
we give an illustration of the two phases that occur in 
BFMs with spinless fermions. 

As shown in Appendix IA1 the scaling exponent of the 
fermionic polaron pairing phase can be calculated to be: 

a f - PP = 2-2[X 2 c K- 1 +K- 1 - 2X c K~ e 1 ] 

= 2-2Kp\ (50) 

where we used the same polarization parameter A c as in 
the previous section, and the same algebraic relation be- 
tween the scaling quantities as for the derivation of (|46[) . 
Note that the expression (|5TJ)) is dual to the scaling of 
the CDW operator, Eq. (|Ig|), and, hence, both the CDW 
phase and the /-polaron phase are governed by a single 
Luttingcr parameter, Kp, which also appears in the sin- 
gle particle correlation function of /-polarons as shown in 
Eq. (gBJ). The fact that the scaling exponents (JUJ), gg} 
and ([50]) are the ones of a Luttinger liquid of ID fermion 
with parameter Kp indicates that the ID BFM system 
should be understood as a Luttinger liquid of polarons. 
We note that fl >x Jn,\ c = fl,\jR,x c : as can be seen 
from the definition of /, Eq. (|4"2")l . Therefore, the scaling 
exponent of Ocdw is not affected by the screening cloud. 
In Fig. DJa) we also show the scaling exponent of the new 
type of order parameter, Of^pp. We note that ctf-pp is 
always larger than apFP, showing that BFP cannot be 
the most dominant quantum phase in the whole parame- 
ter regime of a BFM. Furthermore, the scaling exponents 
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CDW 




FIG. 2: Illustration of the two phases that occur in a BFM 
with spinless fermions, CDW and /-PP. In the CDW phase 
the system develops a 2fcj density modulation in both the 
fermionic and the bosonic liquid. In the /-PP phase, fermionic 
polarons pair up and form a superfluid state. 



shown in Fig. [Ha) demonstrate that divergencies of the 
CDW and /-PP susceptibilities are mutually exclusive 
and cover the entire phase diagram up to the phase sep- 
aration (PS) regime, which can also been seen directly 
from the scaling exponents in Eqs. (J48J) and ([50)) . This 
explains why a single impurity can be irrelavant when 
the system is in the regime between the CDW phase and 
the BFP phase, as discussed above. 

In Fig. [5] we show a schematic representation of a 
CDW and an /-PP phase. In Fig. [IJb) we show a global 
phase diagram of a BFM by considering the FP cou- 
pling (g) and the effective fermion-fermion interaction 
(G) as independent variables. The shading density indi- 
cates the strength of the phonon cloud of the /-polaron 
pairing phase, 2A C . We find that for a fixed effective 
backward scattering between fermions, i.e. G =constant, 
CDW phase, /-PP phase and phase separation regimes 
become dominant successively when the long- wavelength 
fermion-phonon coupling (<?) is increased. 



C. Polaronic effects on bosons 

The polaronic construction, demonstrated and dis- 
cussed in the previous sections for the fermionic atoms, 
can also be done for the bosonic atoms, see Eq. ([^5]) . In 
Fig. [DJa), we also show the calculated scaling exponent 
of the bare boson condensation field, Obb(x) — b(x), 
and that of the bosonic polaron (6-polaron) condensa- 
tion field, Ob-p(x) = b(x). In one dimensional sys- 
tems, the elementary excitations of the fermions are Lut- 
tinger bosons , which can provide the bosonic screen- 
ing clouds around the bosonic atoms, leading to a higher 
scaling exponent than the bare bosons as shown in Fig. 
[lja). However, in the fast boson limit (vj, 3> Vf) that we 
consider here, the dressing of the bosonic atoms is very 
weak (i.e. olbb ~ otb—p) unless the system is close to the 
phase separation region. 




FIG. 3: (a) Phase diagram for a mixture of bosonic and 
spinless fermionic atoms in a ID optical lattice, for the first 
example discussed in Section MI Dl Shading in the /-PP 
phase describes the strength of the bosonic screening cloud 
(2A C , see Eq. (|42p ) around a pair of fermions. Xl and E R 
are the lattice period and recoil energy, respectively. Other 
parameters used for this figure are (see text for notations): 
u b = 4, v f = 0.5, V b ,± = V f ,± = 20E R , V M = 2E R , boson- 
boson scattering length au> = O.OlAz,. (b) Phase diagram 
for a mixture of bosonic and spinless fermionic atoms in a 
ID optical lattice, for the second example in Section MI Dl 
Shading in the /-PP phase describes the absolute value of 
the strength of the bosonic screening cloud (|2A C |) around a 
pair of fermions. Other parameters used for this figure are 
(see text for notations): Uf, = 3, Vf = 0.2, Vf.±/Ef — 20, 
Vb,±/Et — 1.65Vf,x/Ef, boson-boson scattering length abb = 
lOOao, where ao is the Bohr radius. 



D. Phase diagram in terms of experimental 
parameters 

In the last sections we derived and discussed the T = 0- 
phase diagram of BFMs. We will now give two plots 
of the same phase diagram in terms of parameters that 
are experimentally more accessible. In the first exam- 
ple (depicted in Fig. [3] (a)), we assume that the lattice 
strengths Vfrtj/ii can be tuned independently, and we 
keep the densities Vf and Vb, and the boson-boson scatter- 
ing length abb constant. We now vary the boson- fermion 
scattering length abf and the lattice strength V&m of the 
laser that is parallel to the system direction and couples 
to the bosons. This lattice strength mostly affects the 
bosonic tunneling tb'. For small V^y, the tunneling am- 
plitude is large, and therefore the phonons are fast. For 
large V&n the phonons are slow and the system develops 
CDW QLRO. For large values of abf the system experi- 
ences the Wentzel-Bardeen instability, as shown in Fig. 

EDb). 

As a second example we will now discuss a 87 Rb- 40 K 
mixture in an anisotropic optical lattice created by a set 
of Nd:YAG lasers. In Fig. [3](b) we show a phase diagram 
for a mixture of bosons and spinless fermions as a func- 
tion of the scattering length between bosons and fermions 
(abf (in units of ao)) and the strength of the longitudinal 
optical lattice for fermionic atoms (Vf t \\/Ef). We assume 
that the lasers are far detuned from the transition ener- 
gies, and that the optical lattices for the bosons and the 
fermions are created by the same set of lasers. In this set- 
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0.5 1.0 1.5 

g/v f 



FIG. 4: Phase diagram of a BFM with spinless fermions at 
half-filling. All parameters are as in FigfT] b. For large G a 
charge gap (CG) appears due to the Umklapp scattering. 



up, the bosonic and fermionic lattice strengths are now 
constrained to the ratio: Vb,\\/(±)/ Eb = 1.65V/j|/(j_)/.E/ 
for a laser wavelength of 830nm [22| . We chose nega- 
tive values for the boson-fermion scattering length, be- 
cause for most combinations of hyperfine states the scat- 
tering length between 87 Rb and 40 K is negative. Due 
to the different scaling of Vb and Vf with Vb , II /-Eft an d 
Vf^/Ef the ratio of Vb/v/ varies going along the verti- 
cal axis. The lower portion of the diagram corresponds 
to slower bosons (compared to the fermions) the upper 
part to faster bosons. For relatively weak boson-fermion 
interactions and strong confinement the system is again 
in the CDW phase, in which the densities of fermions 
and bosons have a 2fc/-modulation. In the case of very 
strong boson-fermion interactions the system is again un- 
stable to to a Wentzel-Bardeen instability, which now 
leads to collapse (CL) [UHij], rather than phase separa- 
tion. The two regimes are separated by an /-PP phase 
with an increasing value of the polarization parameter A c 
for stronger interaction. 



1.5 



SF 




MI \ 


PS 



0.1 0.2 0.3 0.4 
g/\f 



FIG. 5: Phase diagram of a BFM with unit bosonic filling. 
For small Kb and small g the bosons undergo a Mott insulator 
transition. In this diagram we only indicated the phases of the 
bosons. The fermions will be in an /-PP phase throughout 
the diagram. 



A. Charge-gapped phase for half-filling of the 
fermions 

In this section we still assume the system to be in the 
regime of slow fermions (vf <C Vb) and large bosonic fill- 
ing {yb 2> Vf), but the fermionic filling is set to be at 
half-filling, i.e. Vf = 1/2. In this situation, the derivation 
given in Eq. (fl"6"|) , leading to the backscattering term Eq. 
(f3"2")l . has to include Umklapp scattering, which involves 
a momentum transfer of Akf = Airvf = 2tt. Similar to 
the derivation of Eq. (|32|) . Umklapp scattering can be 
obtained to be: 

In bosonized representation, this expression becomes: 



IV. PHASE DIAGRAMS IN OTHER 
PARAMETER REGIMES 

In the previous two sections we have considered BFMs 
in the limit of large bosonic filling (z/& 3> vf) and fast 
bosons (vb ^> Vf). We also assumed that both the 
fermionic filling and the bosonic filling are incommen- 
surate to each other and to the lattice, so that higher 
harmonics of the bosonization representation could be 
neglected. We will now consider two cases in which these 
assumptions are relaxed in different ways: In section 
IIV A[ we consider a BFM with fermions at half-filling. 
At this particular fermionic filling, phonon-induced Umk- 
lapp scattering can open a charge gap. In Sec. IIV Bl we 
consider the case when the boson filling fraction is unity, 
so that the system can undergo a Mott insulator transi- 
tion. A further discussion of commensurate mixtures is 
given in [32[. 



H uk = -7^2 J dx cos [40/ (a;)]. (52) 

2 

32 k 

(73 is the Umklapp parameter given by 53 = — — , which 

is equal to the backscattering parameter G. Since this 
is a nonlinear term that cannot be diagonalized, we use 
a scaling argument to study its effect. Such a renor- 
malization flow argument at tree-level corresponds to a 
systematic expansion around g% — 0. The scaling dimen- 
sion of this term is 4K^, so the flow equation for is 
given by: 

^ = (2-4^)03 (53) 

Therefore, the Umklapp term is relevant for Kp < 1/2. 
When (73 becomes relevant, the fermion excitations be- 
come gapped due to the strong repulsion mediated by 
the 2fc/-phonons. As indicated in Fig. 01 we find that 
such a charge gap (CG) opens in the large G regime. 
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B. Mott insulator transition for bosons at unit 
filling 

In this section we determine how the Mott insula- 
tor (MI) transition in the bosonic sector is affected by 
the presence of the incommensurate fermionic liquid. 
However, the results of this section equally apply to a 
boson-boson mixture because the dynamical properties of 
fcrmions and bosons in ID are very similar. The starting 
point for our discussion is not the single band Hubbard 
model (Eq. [1]), but a continuous ID system, to which we 
add a weak 2fc;,-periodic external potential: 



2k,, 



= I dxV 2 k b cos26> 6 . 



So the entire system is described by 



H = H, 



(0) 



H 



(0) 



(i) 

bf 



H, 



21,,, 



(54) 



(55) 



For simplicity we neglect the fermion backward scattering 
term. We again use a scaling argument to determine the 
effect of the non- linear term (f5"4"]). This term has a scaling 
dimension of K$, and its flow equation is therefore given 
by: 



2k,, 



dl 



(2 - K s )V 2kb 



(56) 



The condition for the Mott insulator transition is there- 
fore Ks < 2. Since we always have K$ > K},, the presence 
of the other atomic species tends to 'melt' the Mott in- 
sulator. This can also be seen in Fig. [5] For g = we 
obtain a Mott insulator for Kt < 2 and a superfluid for 
Kb > 2. However, when we increase the fermion-phonon 
coupling g the Mott insulator regime shrinks and even 
vanishes entirely before phase separation is reached. This 
is reasonable because the coupling to the second species 
of atoms induces an attractive interaction, which tends 
to decrease the tendency of the bosons to form a Mott 
insulator at unit filling. 



V. SPINFUL SYSTEM 

In this section we consider BFMs that have fermions 
with two internal hyperfine states. We assume the 
fermionic sector to be SU(2) symmetric, that is, the two 
Fermi velocities and their filling fraction are equal, as 
well as their coupling to the bosonic atoms. We again 
consider the limit of fast bosons, vj, ^> Vf. 



A. Hamiltonian 

Analogous to the Hamiltonian of a BFM with spinless 
fermions, Eq. (|33p . we describe the two internal states 
with the fields O-fj and 3>-|-,j.. In terms of these fields, 



the fermionic sector is described by a Hamiltonian of the 
form: 



Bj = 



J i=t,l 



Vf \ Kf ,„ ,9 7T 
7T Kf 



n 2 



+c/ Ti n T nx + u cos(26> T 

11 1 + (2-7ra) 2 v 



26> 



)}, (57) 



where a— >0 + is the intrinsic cut-off of the bosonization 
representation. The parameter U-\{ describes the s-wave 
scattering between two spin states. The last term is 
the backscattering term, which occurs because the two 
fermionic spin states have equal filling. The bosonic sec- 
tor is of the same form as before: 



hS = v 4 



dx 



Kb 

7T 



(58) 



Analogous to Eq. (|3"Tj) . the long wavelength density fluc- 
tuations of the fermions and the bosons are linearly cou- 
pled: 



= u bf J dasn 6 (n t + n|). 



(59) 



Finally, we also integrate out the 2fe/-phonons of the 
bosonic superfluid, leading to the following interaction 
between fermions: 



H ) = E lr / dx [ " 2n * " (5 ^ )2] 

i=T4 77 



(60) 



The total low-energy Hamiltonian is then given by: 



H 



H, 



(i) 

bf 



(61) 



which can be separated into spin and charge sectors 
H = H, 



binations 

e 



p -r Ha, by introducing the following linear com- 
, 1 



p/a = -^=(e T ±e i ), $ P/CT = — (^i^) (62) 



Note that H p coupled to the bosonic field is of the same 
form as a BFM with spinless fermions that we discussed 
earlier and can be diagonalized by a rotation similar to 
Eq. ([35]). 



n b 

n„ = 



= <5iH4 + c$2n Q , 



$ p =7l$A+72$a, (63) 



The diagonalized Hamiltonian of the charge sector has 
the same form as (|3"7|) : 



dx 



j=a,A 



7rU j (x) 2 + -(d x ^ j (x)ri6A) 

7T 
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FIG. 6: Linearized renormalization flow of the parameters 
gi : ± and K a in the vicinity of the SU(2) symmetric fixed point 
at gi : ± = and K a — 1. Systems with SU(2) symmetry are 
forced to be on the separatrix from the upper right quadrant 
to the lower left quadrant. If the initial parameters of such a 
system correspond to a point in the upper right quadrant the 
system flows to the non-interacting, SU(2) symmetric fixed 
point at the origin. For initial parameters in the lower left 
quadrant the system flows towards a strongly interacting fixed 
point and opens a spin gap. 



= (2~2K a ) 9l . ± (69) 

These equations are a systematic expansion around the 
non-interacting fixed point at gi,± — and K a = 1. In 
Fig. [51 we show the linearized flow generated by these 
equations, in the vicinity of this fixed point. For K a > 1, 
the parameters of the system flow towards the fixed point 
(gi.± — 0,K a — 1), whereas for K a < 1, they flow to- 
wards the strong coupling fixed point — > ±oo, K a = 
0). The nature of spin excitations of the ground state fol- 
lows from the well known properties of the sine-Gordon 
model for SU(2) symmetry: When the initial <?u_ is pos- 
itive the system has gapless spin excitations, i.e. gu_^0 
and K a — >1, and when the initial gu_ is negative the sys- 
tem has a spin gap, i.e. gi±—* — oo and K a — >0. 

C. Scaling exponents 



and, similarly, the eigenmode velocities, va and v a , are 
of the same form as 



U a/A 



K 2 



± 



\ \J K - v p) 2 + l ^9 2 p v b v p l^) 



where we defined v p = (v 2 — 4G 2 ) 1 / 2 and g p = g p e e , with 

e 6 = ((v p - 2G p )/(v p + 2G p )) 1 /\ g p by g p = V2g. v p is 
given by v p — Vf + G, G p = G + Gn/2. Gn is given by 
Gii = Uh/2-k. All of these parameters and expressions 
can be immediately obtained by transferring the results 
from the spinless BFM case. 

The quadratic part of the spin part of the Hamiltonian 
can be diagonalized by using 



IT 

where y ' K„ = e 
V a = Vf - Gii, G, 



(66) 



and tanh20,j = —2G a /v cr , with 
= G - Gn/2 and v„ = y/v*- AG% 
The transformed spin sector then becomes a sine-Gordon 
model, 



dx 



nU a (x) 2 + - (d x ^(x)) 



(2na) 



dx cos 



(67) 



Here, gi± = U~n — AirG is the effective backward scatter- 
ing amplitude for fermions. 



B. Renormalization flow 

In this section, we deal with the spin Hamiltonian by 
performing an RG calculation, similar to the calcula- 
tion in section IIVI As discussed in [l(| , and as is well- 
established in the literature, the RG flow at one-loop is 
given by: 



dK a 
dl 



2 



\ TTVn- ' 



(68) 



Now we will determine the scaling exponents of the 
order parameters that can show QLRO in the phase 
diagram. The order parameters of a standard SU(2) 
symmetric LL of fermions are the charge density wave 
(CDW), spin density wave (SDW), singlet and triplet 
pairing (of bare fermions), and Wigner crystal (WC) op- 
erator. When such a system is coupled to a bosonic su- 
perfluid we again need to replace the pairing of bare 
fermions by pairing of polarons (analogous to Section 
LH1 see also Appendix IA")) . First we consider the 2k /- 

mode of the CDW operator (Ocdw — J2 S /J l/s,k) an d 
of the SDW operator (Osdw = £ S)S / /JlO's.s' 
In bosonized form these correspond to Ocdw. ,2k , ~ 
exp[^/2iQ p + \/2iO a ] and OsDW,2k } ~ exp[y/2i<d p + 
\f2i<& cy \. Their scaling exponents are given by 



acDW = 2 - K fj - K a 
olsdw = 2 - K (i - K~ 



(70) 
(71) 



For the polaron pairing phase, both singlet (SPP, 
Ospp = fl,L~f UR - fl,Lfl,n) and tri P let ( TPP > °tpp = 
fl,hfi,R) pairing needs to be considered. In bosonized 
form, these operators are given by Ospp,k=o ~ 
cxp[-2iA c $6 + \/2i$>p + \f2iQ a ] and Otppm=o ~ 
cxp[— 2iA c $h + V2i$ p + V2i$ a ], with scaling exponents: 



a SP P = 2-iK- 1 
a TPP = 2 - [Kp 1 



These expressions were 



by using A c 



(72) 
(73) 



K f K~ 2 , 



obtained 

K e /(V2K je ), and the identity K' 1 
as we did for spinless BFMs. Again, we can see that 
these many-body order parameters have scaling parame- 
ters controlled by two parameters only, Kp and K a , in- 
dicating that the spinful BFM system can also be under- 
stood as a LL of polarons. 

We now discuss how these scaling exponents behave 
in different parameter regimes. We consider the case in 
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FIG. 7: Phase diagrams for a mixture of bosonic and S — 
1/2 fermionic atoms with Vb/vf = 5 and Kt = 10. In (a) 
Ufi/vf = — 0.8w and in (b) Ufi/Vf — 0.8n. Parentheses (. . .) 
indicate subdominant phases. 



which the spin sector flows towards the Gaussian fixed 
point, K a = 1. Therefore, SDW and CDW as well as 
TPP and SPP are degenerate at the level of algebraic 
exponents. However, by taking into account logarithmic 
corrections 10], SDW and TPP tur n out to have slower 
decaying correlation functions and are therefore domi- 
nant. As a result, in the gapless phase there are two 
regimes of QLRO: SDW (CDW) and TPP (SPP), where 
the brackets refer to subdominant ordering. For large G, 
the system is in the spin-gapped phase, because gi t ± is 
negative now, and therefore becomes relevant. The sys- 
tem then flows to the strong coupling fixed point and 
K a — ► 0, according to the RG flow. 



D. Phase diagram 

In Fig. [7] we show two examples for different values 
of f/f,4.. In Fig. (a) we show the phase diagram for 
U hl /v f = -0.8tt, and in Fig. [7] (b) for U hl /v f = 0.8tt. 
As described above, the system develops a spin gap for 
gi± < 0. The phase transition occurs at AttG = £/-f.l- In 
Fig. [7]b), this corresponds to G/v f = 0.2, in Fig. [3(a), 
the entire phase diagram is in the gapped phase, because 
U^i < 0, and G can only be positive. 



E. Wigner crystal 

Apart from these phases the system can also show 
Wigner crystal (WC) ordering. Within the WC phase, 
the fermionic atoms crystallize in an alternating pattern 
of spin-up and spin-down atoms, giving rise to a 4fc/- 
density modulation. The Wigner crystal order parameter 

is given by O wc = J2 S /r, s /.r ,- s /i,-s/i, s ~ exp[y/8iQ p ] 
and has the scaling exponent 

a wc = 2-AK p (74) 

Within the gapless phase (K„ = 1), the Wigner crystal 
phase can become dominant for Kp < 1/3, which can be 
achieved for large values of G. 



VI. EXPERIMENTAL ISSUES 

In this paper we have calculated the phase diagrams 
of different types of BFMs in various parameter regimes. 
In this section we will discuss issues concerning the ex- 
perimental realization of our systems. 

Throughout the paper we concentrated on infinite sys- 
tems at zero temperature, T — 0. The realization of ID 
Luttingcr liquids of ultracold atoms will of course be in 
a finite lattice (around 100 sites) and at finite temper- 
ature. At finite temperature the correlation functions 
become approximately (O(x)O(0)) ~ (£ sinh(a;/£)) Q ~ 2 , 
where £ ~ Vf /ksT is the thermal correlation length. For 
a finite system of length L = N\l (N being the num- 
ber of lattice sites), the T = 0-properties of the system 
are visible for £ ~ L. This corresponds to a temperature 
regime of T ~ Tf/N, with Tf being the Fermi temper- 
ature. So for a system with N = 100, T ~ 0.017/ is 
required. 

Another issue is the experimental signature of the dif- 
ferent phases. Here we present several approaches that 
can be used to detect the quantum phases discussed 
above. To detect the CDW phase, one can perform 
a standard time-of-flight (TOF) measurement. In the 
CDW phase the fcrmion density modulation will induce a 
2k f density wave in the boson field in addition to the zero 
momentum condensation so that the CDW phase can be 
observed as interference peaks at momentum k — 2k t 
in a TOF measurement for bosons, whereas the /-PP 
phase will show a featureless superfluid signature. We 
note that the boundary of an atomic trap and other in- 
homogeneities in a realistic experiment can pin the CDW 
phase and generate a true density modulation. One may 
also use a laser stirring experiment [34| to probe the phase 
boundary between the insulating (pinned by trap poten- 
tial) CDW and the superfluid /-PP phase: a laser beam 
is focused at the center of the cloud to create a local po- 
tential, and is then moved oscillatory in the condensate. 
If the system is in the pairing phase, the laser beam can 
be moved through the system without dissipation, which 
would manifest itself as heating, if only its velocity is 
slower than some critical value [34| . At the /-PP / CDW 
phase boundary this critical velocity goes to zero, re- 
flecting a transition to the insulating (CDW) state. This 
scenario follows from the RG analysis of a single impurity 
potential [37j], as described in Appendix [Cl To probe the 
boundary of phase separation one can measure the dipo- 
lar collective oscillations of the system, generated by a 
sudden displacement of the harmonic trap potential with 
respect to the lattice potential [36]. When the system 
is near the PS boundary, fermion-boson interaction will 
reduce the frequency of the dipolar mode essentially to 
zero, because the two atomic species become immiscible. 

As the most promising approach we propose to study 
the noise correlations of TOF measurements, as discussed 
in 35]. This type of measurement treats particle-particle 
correlations (i.e. pairing fluctuations) and particle-hole 
correlations (i.e. CDW and SDW fluctuations) on equal 
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footing, and therefore reflects the formal duality of these 
two types of phases accurately. Furthermore, this ap- 
proach seems to be well-suited for the study of Luttinger 
liquids, because it gives signatures of the various fluctu- 
ations in the system, that dominate in different regimes 
of the phase diagram. 

VII. CONCLUSION 

In summary, we used bosonization to investigate the 
quantum phase diagrams of ID BFMs. Interactions be- 
tween atoms can lead to interesting phenomena like spin 
and charge density waves, singlet and triplet pairing of 
atoms. We introduced polarons, i.e. atoms of one species 
surrounded by screening clouds of other species, and ar- 
gued that the rich phase diagrams of BFMs can be nat- 
urally interpreted as Luttinger liquid phase diagrams of 
such polarons. We also considered several commensurate 
filling cases and obtained gapped phases in some parame- 
ter regimes. We discussed several techniques for probing 
our results experimentally. 

We thank E. Demler, M. Lukin, H. Moritz and B.I. 
Halperin for useful discussions. 

APPENDIX A: SCALING EXPONENTS 

In this section we will derive the scaling exponents of 
a number of operators that were considered in the search 
for the phase diagram. In section [AT] we discuss BFMs 
with spinless fermions, and in section IA 21 BFMs with 
spinful fermions. 

1. Spinless Fermions 

For the case of s pinl ess fermions, a broad class of op- 
erators 0(x) (cp. Jl6|) can be written in terms of the 
Luttinger fields as 

0{x) ~ exp[ Y. ('A,. : B, • 'A, i (Al) 

j=bj 

or a sum of products of this type. ^j,i/2 are arbitrary 
real parameters. To derive the correlation function of 
such an operator we first write 0&// and $6// in terms of 
the eigenfields &A/a and &A/a, as given by the transfor- 
mation (|35|) . The correlation function of the above given 
general operator, Co(x) = (O(x)O(O)) , then behaves for 
large distances as : 

C (x) ~ \x\- a cos((X fA v f + X bA u b )wx), (A2) 
with the scaling exponent a given by 

+AL^r' + ^h*- 1 + 2A M A / , 2 if- 1 ) (A3) 



From this general expression we can deduce the scal- 
ing exponents discussed in this paper, by specializing the 
operator 0{x) to various cases. 



2. Spinful fermions 

For the case of spinful fermions, we consider an oper- 
ator 0(x) given by: 

O(x) ~ exp[ Y : 'A,,. .<-)., • 'A ; j<!>, (A4) 

j=b,p,a 

To determine the correlation function of this operator, we 
proceed as in the previous section. First, we use the linear 
transformation (|63p . to write 9 p / b in terms of the fields 
^ A/ a- Since the density/boson-sector of the Hamiltonian 
is just a sum of Gaussian models in terms of these fields, 
the scaling exponents can be determined as in IA 11 For 
the spin sector we use the RG results discussed in Section 
IV Bl If the non-linear term ir irrelevant the system is 
described by a Gaussian model with K a = 1 . If the non- 
linear term is relevant we have K a — > 0, i.e. any overlap 
of the operator with the phase field $ CT will make the 
correlation function short-ranged, whereas the field 9 a 
acquires ordering. 

The correlation function of the operator O(x), 
Co{x) = (O(x)O(0)), can be determined to be: 

C (x) = \x\~~ a cos((A/,i2// + \ byl v b )Trx) (A5) 

The scaling exponent a is given by 

a = -(\l tl Ks + >? p ,iKp + 2A{, ! iA P: i if^ 

+ A^tf- 1 ) (A6) 

In Table [TTl we first give the scaling of the single-particle 
operators b and f-\/i, as well as the standard order pa- 
rameters of a LL of spinful fermions, i.e. fjf s and /j/j. 
(the 2k f -modes of these operators correspond to CDW 
and SDW ordering, respectively), the pairing operators 
/t/j. (singlet) and f s f s (triplet), and the order parameter 
of the Wigner crystal, /|/|/x/|- Apart from these stan- 
dard order parameters of LL theory, we consider a wide 
class of operators of the form (/ T f ) m /f {ft) p f?(tf) r b s . 



APPENDIX B: POLARON EFFECTS 

In this subsection we will verify that the conventional 
construction of polaron operators based on the canonical 
polaron transformation (CPT) [3l[ is equivalent to the 
construction within the bosonization approach presented 
in Section IIII Al of this paper. We will also discuss the 
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operator 0(x) 


wavevector q 


exponent a 


remarks 


b 


2mfc;, 


i 

2 


'(2m) 2 K s + K- r 




f 


(2n + l)k f 






Pf 


2nkf 


2n 2 K f3 


n / 


ff 


Inkf 


2n 2 K p + 2K- 1 




tfb 


2mk), 


2m 2 K s 


m / 


[ & t]mj 


2mkt > 


l 

2 


~{2m) 2 K s + mfKr 1 ' 








2nkf 


1 

2 


\2n) 2 K () + (2n 1 ) 2 K~ 1 ' 






2mkt, + 2nkf 


1 

2 


'(2m) 2 Kg + (2n) 2 K f3 - 2(2m)(2n)K f3S 


m,n/0 




2mkt + 2nkf 


i \(2m) 2 K 5 + (2n) 2 K fi - 2(2m)(2n)K l35 + m\K- r 




[&t] mi 


2mkt + 2nkf 


| [(2m) 2 Kg + (2n) 2 K - 2(2m)(2n)K 0S 
+TTIIK- 1 + (2n 1 ) 2 K~ 1 - 2m 1 (2n 1 )K- 1 ] 




[fot]" 11 [/] 2 ™1 +1 


2mk b + (2n + l)kf 


| [(2m) 2 K s + (2n + l) 2 K fj - 2(2m)(2n + l)K 0e 
+m\KT 1 + (2m + lfK' 1 - 2mi(2m + l)^ 1 ] 




/ 


(2n + l)k f 


2-\\{2n + l) 2 K f) + Kg 1 } 




ff 


2nkf 


2n 2 Kf3 + 2K^ 




b 


2mkt > 


2-| \(2m) 2 K s + K^ 







TABLE I: Table of Luttinger exponents and characteris- 
tic wavevectors of various operators composed of boson and 
fermion operators, m and n are arbitrary integer number re- 
sulting from the higher order harmonics of the bosonization 
representations. 



behavior of the polarization parameter A c . The CPT and A specifies the strength of the phonon 

operator is given by dressing. 

Ux = e -*£^o( F ^"/,*- h - c -) (Bl) 

where (3k is the phonon annihilation operator and n/^ 
is the fermion density operator. Ff. is given by Fk — 



The sum over the wavevector k is over the regime of 
acoustic modes. When applied to a fermion operator, 
the CPT gives 

u^f(x)u x = /Wf ; A EM.(^^^ h -> ( B2) 

which is the standard expression for a fermionic polaron 
operator [3l[ . In terms of Luttinger bosons the fermion 
density nf t k is given by: 

n f , k = ^/kL/2^(B k +Bl k ), (B3) 
and therefore U\ becomes: 

Ux = ^E^MMt^^) (B4) 

When applied to a fermion operator, this transformation 
gives 

U^f(x)U x = f{x)e~ iX ^ (B5) 



because U^^fUx = — A$f>, which can be seen from 
the operator representation and ([50]) . and by apply- 
ing the Baker-Hausdorff formula. The superficial differ- 
ence between Eq. (|B2[) and Eq. jB5| can be overcome by 
a trivial phaseshift (3k — * isgn(fc)/3fc. As discussed in this 
paper, A c = K e /K je is the most appropiate choice for 
the fermionic polarons, because the correlation function 
of such polarons has the slowest algebraic decay. In Fig. 
|S] we plot this polarization parameter A c in comparison 
to the value of 'complete dressing' Ao: 

A = ^/Kfc (B6) 

Vb 

To obtain this quantity, which is the 'size' the polariza- 
tion would have if the fermions were static, we set the 
Fermi velocity to zero (vf — 0), so that the Hamiltonian 
of the system becomes: 

H = J2 v b \k\Plh + £ g\k\(/3{ + P- k )(B k + BIJB7) 

k^O k^O 
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operator 0(x) 


wavevector q 


exponent a 


remarks 


b 


2mfcj, 


^mfKg + K- 1 ] 




fr/i 


(2n + 1)% 


(n + i) 2 (if/3 + if,) + + if" 1 ) 




fif* 


2nfe/ 


n 2 if,3 + n 2 if CT ] 


n / 




2(n + m)kf 


\n + m) 2 ^ + (n - m - \) 2 K a + if" 1 ] 


n / 


fsfs 


2nkf 


n 2 Kp + n 2 K a + if" 1 + if" 1 " 








2(n + l + l)k f 


\n + l + 1) 2 K P + K- 1 + (n - Z) 2 if CT ~ 






fjflhfi 


2(n + m)kf 


"(n + m) 2 Kp + (n - m) 2 Ko\ 






2mkt > 


2m 2 K s 


m / 




2mfc;, 


'2m 2 K s + \mlK- 1 } 




lfs] 2ni 


2nkf 


n 2 (Kp + K a ) + njiK- 1 + K' 1 )] 




[/s] 2„ 1+1 


(2n + l)fc/ 


\n + D 2 {Kf3 + K a ) + (m + + if- 1 )] 






2mfc;, + 2nfc/ 


[2m 2 K s + n 2 K(3 - 2-j2mnK ps + n 2 K a 
+ \m\K7 l + nlK' 1 - s/2m x n x K-} + nfif" 1 ] 




[fo f ] ml [/] 2nI+1 


2mfc 6 + (2n + 1)*;/ 


[2m 2 if 5 + (n + \fK - V2m(2n + 1)K 0S + (n + ±) 2 if CT 
+ ^m 2 1 KZ 1 + (ni + IfK- 1 - V2m x (n x + ^if" 1 + (m + i) 2 *^ 1 ] 




[&T 1 [/t/x] 


2mfci, + 2(n + / + 


[2m 2 if a + (n + Z + if Hp - 2^/2m(n + I + l)K s + (n - Z) 2 if CT 
+ im 2 if e " 1 + K- 1 - ySmiif" 1 ] 




b 


2mfci, 


i[(2m) 2 if a + if,- 1 ] 




fr/i 


(2n + l)fc/ 


(n + i) 2 (if,3 + if CT ) + f(ifg 1 + if" 1 ) 




fjs 


2nkf 


n 2 if,3 + n 2 if CT + if- 1 + if" 1 " 






hh 


2(n + l + l)k f 


"(n + Z + 1) 2 K + Kp 1 + {n- l) 2 K a ~ 







TABLE II: Table of Luttinger exponents and characteristic 
wavevectors of various order parameters composed by boson 
and fermion operators for spinful system, m, n and Z are 
arbitrary integer number resulting from the higher order har- 
monics of the bosonization representations. 



In this limit of infinitely heavy fermions, a term of the 
type (fT?)) does not exist. In (|B7[) the linear term simply 
shifts the bosonic modes by an amount g/vb(Bh + B_ k ). 
As we can see from (|B4[) this corresponds exactly to the 
polarization Ao, Eq. (|B6[) . However, for finite hopping, 
or finite mass of the fermions, this value is reduced, be- 
cause the polarization cloud cannot entirely 'follow' the 
fcrmionic atoms. For Ubf — * one finds 

A c /A -> — H*_ (B8) 

Vb + Vf 

In Fig. [8] we show A c plotted for different sets of param- 
eters. It smoothly interpolates between the small inter- 
action limit (|B8jl and the static limit (|B6j) . The static 
limit is achieved in the vicinity of phase separation be- 
cause the fermionic effective mode velocity Va goes to 

I 



APPENDIX C: RG FLOW OF AN IMPURITY 
TERM 

We consider the following local impurity term in our 
Hamiltonian: 

H imp = [ dxV{x)f{x)f{x) (CI) 



zero when the system approaches separation. This corre- 
sponds to infinitely heavy fermions, and therefore to the 
static limit. 




FIG. 8: Plot of the 'dressing' parameter A c compared to com- 
plete dressing Ao, for iff, = 5 and G/vf = 0.1 as a function of 
g/v f . With v b /v f = 2, 3, 5 for I-III. 



I 

with V(x) strongly spatially peaked around at x = 0. 
This term leads to an additional term in the action given 



16 



by 

S imp = Y, V -Y f dTe M^mO f ) (C2) 

m 

where v m is the Fourier transform of the potential V(x) 
around 2mkt. The RG flow of these terms at one- loop, 
as discussed in [37[ and [33| is given by: 

^ = (l-m 2 Kp)v m (C3) 



Therefore, such an impurity term is relevant exactly for 
Kp < 1, and irrelevant outside of the CDW phase. The 
physical interpretation of this result is as follows: In the 
CDW regime a local impurity 'pins' the charge ordering 
of the system. Outside of the CDW phase, the irrelevance 
of an impurity term in the system is an indication of 
a superfluid phase, which must be provided by polaron 
pairing. 
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